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Abstract. Nuclear electron capture rate from continuum in an astrophysical plasma 
environment (like solar core) is calculated using a modified Debye-Hiickel screening 
potential and the related non-Gaussian g-distribution of electron momenta. For q = 1 
the well-known Debye-Hiickel results are recovered. The value of q can be derived from 
the fluctuation of number of particles and temperature inside the Debye sphere. For 
^Be continuum electron capture in solar core, we find an increase of 7 - 10 percent 
over the rate calculated with standard Debye-Hiickel potential. The consequence of 
this results is a reduction of the same percentage of the SSM *B solar neutrino flux, 
leaving unchanged the SSM ''Be flux. 
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1. Introduction 

Since the early works by Bethe [1] and Bahcall [2], [3] great attention has been devoted 
to the screening effect of Coulomb potential on electron capture (EC) by nuclei in 
astrophysical plasmas [4] and on its implications with neutrino production. 
The plasma influence has been in great part explored by means of Debye-Hiickel (DH) 
treatment. However, Johnson et al. [5] found that assumptions for the validity of 
DH potential are strongly violated in stellar cores, because, among other reasons, the 
requirement to have many particles in a Debye sphere is not fulfilled. At solar conditions 
there are only few (about four) particles per Debye sphere. Role and importance of 
several effects in the electron component of a weakly nonideal hydrogen plasma have 
been recently investigated in [6], inducing us to investigate on EC beyond DH approach. 
The contribution to nuclear EC rate of bound and continuum electrons has been widely 
studied by Shaviv and Shaviv [7J. They have found that under solar conditions all 
elements with 4 < Z < 12 are fully ionized, therefore the contribution from bound 
electrons to the capture rate can be disregarded. Formerly, Bahcall and Moeller [8] 
reported that the ratio of the rate from continuum over the rate from bound states is 
about 1.2. 

Fluctuations of electric microfields and Debye particle numbers have been investigated 
by Gruzinov and Bahcall [9] and by Brown and Sawyer [10] within the Feynman density 
matrix treatment and the mean field theory applied in the Boltzmann limit validity, i.e. 
at global thermodynamical equilibrium. They have found that non-spherical fluctuations 
can change the reaction rate of about 1-2 percent and that quantum corrections to 
Maxwell-Boltzmann (MB) limit for Z = 4 are less or about 1 percent. 
Screening of Coulomb field of the nucleus by outer electrons has been explored by 
considering the Hulthen potential that has a shape very close to a screened Coulomb 
potential. By solving analytically Klein-Gordon and/or Schroedinger equations one can 
calculate the electron density at the nucleus and consequently derive the Fermi factors 
in terms of hypergeometrical functions [H] . 

In all above quoted papers, continuum EC rates are evaluated in the classical statistics 
limit, using the MB electron momentum distribution that is the correct distribution for 
Coulomb and DH potential V^^{r), but is not the appropriate distribution to be used 
with the modified DH (MDH) potential Vq{r) (\/g=i(r) = V^^{r)) introduced by us in 
[T2] . The spatial charge distribution related to the last mentioned screening potential 
differs from the others having a spatial cutoff. The linear Poisson equation used to 
deduce the screening DH potential must be substituted by a non-linear equation to 
take into account particle correlations and fluctuations and Vq{r) is deduced as a power 
law function. In |12] we have also derived the potential V^(r) in the framework of the 
super-statistics approach [I3l [E], by considering an inverse Debye radius fluctuating 
around its average value given by ^ = (l/i?DH)^- The evaluation of the rate requires 
non-Gaussian generalised g-distribution of electron momenta with the same value of q 
of the spatial charge density distribution. 
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In this work, we evaluate the variation of the EC rate due to the MDH screening 
potential over the pure Coulomb screening rate as a function of average Debye radius 
^ and entropic non-extensive parameter q [151 HBl [E] • We consider the case < g < 1 
to have a screening potential with finite spatial range. This peculiarity of some high 
density astrophysical plasmas and of most of the Sun is due to the fact that, in these 
systems, the mean interparticle distance is smaller or slightly smaller than DH radius 
[7j. Standard DH results are obtained in the limit g 1 of the electron non-Gaussian 
distribution. 

The value of the parameter q can be fixed by means of its relation with the fluctuation of 
particles number and of temperature inside the Debye sphere, therefore the same value 
of q can be used both for spatial charge density and for momentum distributions. 
When this treatment is applied to EC by ^Be in the solar core, we can derive the 
amount of the rate enhancement for the electrons of continuum over the DH rate. 
We find that, at solar conditions, the rate can increase of about 7-10 percent over 
DH screening rate with q ranging between 0.84 and 0.88 (small deformations of MB 
distribution). This result can be useful in the interpretation of the observed ^Be and ^B 
neutrino fluxes, in the evaluation of relevant astrophysical S factors and of CNO solar 
neutrinos [TH [191 1201 [21] . For recent discussions and comments on the EC by ^Be from 
a continuum three-body initial state we send to [22] and references therein. See |23j for 
EC rate evaluated with a percentage of non-thermal fat-tail electron distribution. 
In Section 2, we report the expressions of the rates evaluated with g-generalized 
distribution at few different values of q. In Section 3, we present numerical results of 
the rates, for electrons belonging to the continuum for different values of the parameter 
q and focusing to the case of EC by ^Be in the solar core. Finally, in Section 4, we 
outline our conclusions. 

2. Rates 

It is known that in nuclear continuum EC the evaluation of the rates includes the Fermi 
factor, i.e. the electronic density at the limit of r — > in a pure Coulomb potential. 
Therefore, the pure nuclear rate is corrected because of the Coulomb interaction between 
the captured electron and the nucleus. Usually, within a plasma environment, DH 
potential is adopted in place of Coulomb potential. In this case the electron density at 
the nucleus is known only numerically from the solution of the appropriate Schroedinger 
equation. 

In an astrophysical environment, the momentum distribution of screening electrons 
can differ significatively from MB distribution. We have recently introduced a MDH 
screening potential to take into account non-linear and correlation effects by means of 
fluctuation of 1/^ = {l/R-on) [12]. Of course, also in this case we need to evaluate the 
electron density at r = by solving the appropriate Schrodinger problem for electrons 
in the continuum. 

Whereas Schroedinger solution with a Hulthen potential can be given in a close analytic 
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form, the MDH potential (as well as the standard DH potential) admits only numerical 
solutions for the electron density at r = 0. Alternatively, one can use the Hulthen 
potential that fits quite well the MDH potential in the small r region (near the nucleus), 
but contains an infinite tail instead of having a cut-off at an appropriate value of r = r^nt 
as for the MDH potential with g < 1. A Coulomb cut-off potential, obtained by imposing 
the condition V^{r) = for r > r^ut and consequent discontinuity in the potential at 
= '"cut; was used in the past as the simplest way to screen the Coulomb field |24j . 
Otherwise, in our approach, the cut-off condition arises naturally and the MDH potential 
is a smooth function for any r > 0, vanishing for r > Tcut- 

In the following, we evaluate the rate for the free electron capture by a {A, Z) 
nucleus, given by the integral, in the three dimensional space of velocities, of electron 
capture cross-section (Xe times the electron velocity v, the normalised probability density 
{F^, Fpjj, or Fq) that an electron of the continuum spectrum, with velocity v and 
travelling in a screening potential (Vc, Vdh, Vh or Vg), be at the nucleus with coordinate 
r = and the normalised probability that the electron velocity be v, probability given 
by the distribution function fq{v) {q = I for C, DH and H, where /^^^(f) = /mbI"^) is 
the normalized MB distribution of electrons). This distribution must be appropriate to 
the screening potential used. 

We define the pure Coulomb nuclear electron capture rate, averaged over a MB 
distribution, as 

oo 

7^, (T) = J (ae v) {v) 4 vr dv, (1) 



where 

^e = -^,-(Wo + Wyx, (2) 

is the nuclear electron capture cross section [21 [3] with G the Fermi constant, Wq the 
nuclear energy release for one electron with total energy W, x = Cy (1)^ + Cj^ (cr)'^ the 
well-known reduced nuclear matrix element |2j. 

The Fermi factor for Coulomb potential, given by F^{E) = 2 vr 17/(1 — e~'^'^^) with 
T] = 4/(ao p) where oq is the Bohr radius and p = rUe v is the electron momentum, follows 
from the definition F^{E) = limr^Q\ip(^{r)/pr\'^, where ip^{r) is the wave function of 
the Schrodinger equation with the Coulomb potential. 

The Hulthen rate TZ^ has been evaluated by averaging over MB distribution and 
by substituting in the Coulomb rate of Eq. ([1]) the Fermi factor F^{E) = 
lim^^o |V'h('")/p'"P, where ip^ir) is now the wave function of the Schrodinger equation 
with the Hulthen potential Vk(r) = -(ZcV^dh) e'^''/^^^ (1 - e-^'^/^^^y^ and, in this 
case, the Fermi factor can be obtained analytically [HI [25] . 

Finally, the non extensive rate TZg, can be obtained by substituting in Eq. ([1]) the factor 
F^{E) with the Fermi factor Fq{E) = lim^^o |'^(j('")/P'"P where ipq{r) can be obtained 
as a numerical solution of the Schroedinger equation with the modified DH potential 
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n(r) = -^[l-(l-g)5,r]^/^-^ , (3) 
with 5, = l/[{2-q) i]. 

Consistently with the derivation of the MDH potential in the definition of TZq we must 
insert in place of the MB distribution the normalised non extensive distribution fq{v) 
[12] defined in 

with expg(x) = [1 + (1 — q)x^_1^^'^ the g-exponential {[x\j^ = x6{x) where 6{x) is the 
Heaviside step function) and for q < 1 



5-3g3-gr(i + ^ 



For q = 1 formula (jl]) reduces to the MB distribution. 

The integral for the rate TZq, when g < 1, is performed over the real interval [0, Vcut] 



with fcut = \j2kT/\{\ — g) me] < 1 which defines a cut-off condition in the velocity 
space. When q 1 the rate TZq reduces to DH rate TZj^^. 

3. Results 

We pose our attention to the case of solar EC by ^Be. Results can be easily extended 

to any {A, Z) nucleus (calculations are in progress). 

In figure 1, we report, for several value of ^ the quantity 

where X = H, DH and q (with q = 0.95, 0.85, 0.75) that represents the percentage 
variation of probability density at r = compared to the probability density when the 
screening potential is a pure Coulomb potential. 

Bahcall and Moeller [8] and Gruzinov and Bahcall [9j have found that, for R^^ > 0.4 
and Z = 4, the quantity AP^^{E) is less than 1 percent, for 0.3 Oq < -R^h — 0.4 oq 
about 2 percent and have concluded that plasma screening is unimportant for capture 
from continuum. Our calculation of DH electron density at r = (which corresponds 
to the case of g = 1) agrees with their results. However, in [8] and [9] the above 
authors have calculated the rate averaging over MB distribution, in the frame of a 
global thermodynamical equilibrium, therefore they have found a negligible screening 
effect over Coulomb rate and a negligible effect of fluctuations and correlations over DH 
rate. Although electron density at r = due to Vq{r) is smaller than Coulomb density, 
in the velocity space the probability density in the low momentum region is greater than 
MB because the continuum electron distribution fq{v) we use in this work privileges low 
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Figure 1. Relative variation APy.{E) of the Fermi factor F^{E) over the Coulomb 
Fermi factor {E) for ^ = 0.45 ao and ^ = 0.90 ao and energy keV <E <25 keV. 



momentum electrons. Therefore, screening may be important in continuum EC rate. 
This can be seen in Table 1, where we report the calculated deviations of the rate TZy.{T) 
respect to 7?.^ (T) , by means of the function 

at the value kT — 1.27 keV (where EC by ^Be takes place), for the three values of q 
of figure 1. 
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Table 1. Table 1. Relative variation AVl^{T) of the rate TZ^{T) over the Coulomb 
rate TZ^{T) for several values of ^, at the temperature kT — 1.27 keV, for the three 
values of q of figure 1 . 







H 


DH 


q = 0.95 


q = 0.85 


q = 0.75 




= 0.30 ao 


2.24% 


0.58% 


-2.92% 


-9.45% 


-15.44% 




= 0.45 ao 


1.82% 


0.69% 


-2.65% 


-9.03% 


-15.06% 




= 0.60 ao 


1.62% 


0.45% 


-2.91% 


-9.36% 


-15.48% 




= 0.75 ao 


1.52% 


0.23% 


-3.13% 


-9.58% 


-15.70% 




= 0.90 ao 


1.45% 


0.08% 


-3.28% 


-9.73% 


-15.82% 



For any g < 1, TZq > TZ^. We have verified that deviations depend very smoothly 
on kT except for ^ < 0.45 oq and depend very strongly on q. 

The value of q for EC by ''Be in solar plasma can be derived from the expression that 
links q to fluctuation of [12j. By using the equation of state for g-nonextensive 

systems p6] we obtain for < g < 1 

AA^^„ 1 1 



1 - 9 = -^r^ = I . (8) 

^^oH Jn^ ,/1 + il-q)K 



where K=lN^^ In [0.211 • (0.45 ao)^] - H^j^^-j 



In the solar core, where the average electron density is = 9.1 , ^dh' number 
of particles inside the Debye sphere, is about 4, we can derive q = 0.86. It is more 
safe to consider a range of values of q between 0.84 and 0.88. At kT = 1.27 keV and 
^ = 0.45 ao the calculated TZqiT) is estimated to be about 7-10 percent over standard 
DH (g = 1) estimate. 

Let us consider the ^Be - p fusion, reaction producing '^B and, as a consequence, ^B 
neutrinos, in competition with ^Be electron capture [27] . We have verified that the effect 
of the MDH potential over its rate is negligible. In fact, correction to is effective 
only at relative '^Be - p energies lower than 2.4 keV where fusion cross section has a 
negligible value because its most effective energy is at 18 keV. Therefore, if EC rate 
of ^Be increases over its standard evaluation of a given percentage, ^Be increases its 
destruction while the neutrino flux from ^Be does not change because the ''Be density 
decreases. However, the ^B flux should diminish of the same percentage. This behaviour 
is in line with what is found in experiments |21j . 

4. Conclusions 

We have calculated the EC rate lZq{T) for electrons of the continuum spectrum, with the 
use of the MDH screening potential Vq{r) (derived in p^). The value of the parameter 
g has been estimated from its relation with the fluctuation of the number of particles 
contained in the Debye sphere, in the solar core. It is reasonable to take the range of 
values 0.84 < g < 0.88. We have calculated the rate averaging over the electron non- 
Gaussian g-distribution fqiy). 
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Considering the EC by ^Be at ^ = 0.45 and kT = 1.27 keV we have evaluated an 
increase of the capture rate TZqiT) of 7 - 10 percent over standard DH rate 71^^{T) 
that is the 0.69 percent smaller, at the same conditions, than Coulomb rate TZ^{T). Of 
course, a smaller value of q should imply a much greater enhancement of EC rate over 
DH one. The main consequence concerns the calculated neutrino flux that decreases 
of the same percentage respect to its evaluation with standard DH screening, while ^Be 
flux remains unchanged. 
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